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We consider the supercurrent flow through gated 
mesoscopic semiconductor hetrostructures in which a two- 
dimensional normal constriction is confined between super- 
conducting electrodes. We show that for these structures the 
Josephson current, carried by quantized electron modes, can 
be strongly affected by an electromagnetic field. Photon- 
assisted Landau-Zener transitions between Andreev bound 
states in the constriction manifest themselves in one of two 
ways: (i) If the phase difference between the superconducting 
elements is fixed a series of reversals in the direction of su- 
percurrent flow occurs (Somersault effect), (ii) If instead the 
junction is current biased a series of voltage spikes is seen. 
Both manifestations are a direct consequence of the discrete 
energy spectrum of the mesoscopic junction. We discuss nec- 
essary conditions for the described phenomena to be experi- 
mentally observable. 
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The Josephson coupling between two superconduc- 
tors is associated with charge transfer through the non- 
superconducting region that separates them. The cou- 
pling is therefore determined in a crucial way by the 
nature of the electron states in this region. Quan- 
tum mechanical tunneling through an insulating barrier 
(SIS-junction) §, itinerant propagation of free electrons 
through a normal region (SNS) and propagation of 
plasma waves (boson modes) through a ID-channel of 
strongly correlated electrons || are examples of charge 
transfer mechanisms in different types of weak links. 

The electron states carrying the supercurrent in the 
non-superconducting region can be influenced by an ex- 
ternal time-dependent field. This possibility raises inter- 
esting questions on the nature of Josephson coupling due 
to nonequilibrium electron states. This is the problem 
addressed here. 

Nonequilibrium effects can be expected to be particu- 
larity important in situations where the Josephson cou- 
pling is mediated by a normally conducting microcon- 
striction, where only a few electron states — Andreev 
bound states — carry the current. Weak links through 
ballistic microconstrictions in the two-dimensional (2D) 
electron gas of gated semiconductor heterostructures 
have recently been observed experimentally [Q. An im- 
portant feature of this structure is that the electron en- 
ergy levels in the constriction can be easily influenced by 
applying a potential to the gate electrodes. In this way 



the electron concentration (and hence the Fermi wavevec- 
tor, &p) in the channel can be controlled. The potential 
due to induced charges on the gate electrodes further- 
more provides a mechanism for coupling the 2D electrons 
to an external electromagnetic field. 

The geometry of the system to be considered in the 
following is shown in Fig. 1. The Josephson coupling 
between two superconductors can be expressed in terms 
of the phase difference (f> between their respective order 
parameters. We will first consider this phase difference 
to be a fixed quantitity, which in a SQUID geometry can 
be controlled entirely by the magnetic flux (see inset of 
Fig. 1). Andreev reflection at the boundaries between 
superconducting and normal regions in a SNS configura- 
tion is known || to lead to a discrete set of energy levels 
within the energy gap of the superconductors. Due to 
spatial quantization in the transverse y-direction, the en- 
ergy spectrum of the Andreev bound states in a narrow 
and short constriction consists of a discrete set of pairs 
of levels labelled by the quantum number n M , 



E n>± = ±A Jl - T n {k F ) sin%/2). 



(1) 



Energy is measured from the Fermi level and the trans- 
mission eigenvalue T n (kp) is related to the propagation of 
normal Fermi level electrons through the junction (along 
the x-direction) . For a short junction all of the Joseph- 
son current is carried by these Andreev states. It is 
important to note that states below [E n ,—) and above 
(E n _ + ) the Fermi level carry current in different direc- 
tions. In Fig. 2a, the position of a pair of levels E n> ± 
is indicated by the direction of the corresponding partial 
currents marked by arrows. The dashed arrow for the 
state E n + above the Fermi energy illustrates that this 
state is unoccupied at zero temperature and hence does 
not contribute to the total current. 

We shall consider a situation where a gate-induced 
time-dependent potential is a sum of one part Vb(t), 
which varies slowly on a scale related to the interlevel 
spacing, and another rapidly oscillating part V^cos(o;t), 
where u> is of the order of the interlevel spacing. The 
high-frequency component of the induced potential pro- 
duces interlevel transitions, which result in a strong cou- 
pling of pairs of Andreev states at resonance. Due to 
the change in the interlevel distances resulting from a 
slow time variation of k F the condition for resonance, 
u> = Lu a fj = (E a — Efj)/h (a = n, ±), can be met for any 
pair of levels at some time t = t Ies even if they are out of 
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resonance with the electromagnetic field at t = 0, when 
the potential is switched on. Such a drift into and out of 
resonance is illustrated in Fig. 2, where panel (a) shows 
a non-resonant configuration. When the two levels drift 
into resonance, as in Fig. 2b, resonant interlevel mixing 
results in a finite population of the upper level and a cor- 
responding contribution to the partial Josephson current 
in the reverse direction. 

If the Andreev levels pass slowly enough (see below) 
through the resonance, the result of the dynamic evolu- 
tion is a complete depopulation of the lower level. Hence, 
as the levels drift out of resonance the Josephson current 
has been turned around and the current is flowing in the 
opposite direction (Fig. 2c). After some time out of reso- 
nance, the system relaxes back to its ground state and the 
Josephson current returns to the forward direction. The 
microwave-induced "Somersault" has been completed. 

To discuss the above effects quantitatively we begin 
with the time dependent Bogoliubov-de Gennes (BdG) 
equation: 



im(r,t) = [H + V(t)) ¥(r,t). 



Here H is the standard Hamiltonian for both the nor- 
mal electrons in the microconstriction part of the device, 
\x\ < L/2, (see Fig.l) and for the electrons in the su- 
perconducting regions, \x\ > L/2. The length of the 
junction L is assumed to be small in comparison with 
the coherence length £sm — vf/A (vp is the Fermi ve- 
locity in the constriction). The normal electron Hamil- 
tonian includes an electrostatic potential that confines 
the electrons in the transverse direction and also de- 
scribes the normal scattering from impurities and from 
the barriers at the Sm-S boundaries. The spatial varia- 
tion of the superconducting order parameter due to the 
proximity effect is small and can be neglected since the 
width of the constriction is much smaller than the su- 
perconductor coherence length £s 0. Hence we use 
a step- function model for the pair potential, A(r) = 
A0(|x| — L/2)exp(« sign(a;)<^/2). The time dependent 
potential, V(t) = [V Q (t)+V w cos uit]a z Q (L/2- \x\), where 
o z is a Pauli spin matrix, is confined to the normal region 
and does not mix the electron- and hole-like components 
of the wave function 1 3/. If its high frequency part V u 
is absent, the only effect of the external time-dependent 
field V(t) is an adiabatic energy shift of the Andreev lev- 
els caused by the parametric dependence of E a on kp — 
[{2m/h 2 )(fi- Vb(t)] 1/2 (seeEq. ([[])). If the field has a fi- 
nite high frequency component, different Andreev states 
are mixed due to interlevel transitions described by the 
matrix element V a/3 = V u (^ a (r), a z Q(L/2 - |x|)*^(r)). 
In what follows we will consider the weak coupling limit 
V a /3 <C hu> a p, in which almost all Andreev states are 
only slightly renormalized by the electromagnetic field. 
The exception is when the condition for resonance be- 
tween two states, ui = LJ a p(t), is fullfilled at some time 



t = t TCS . At resonance the interaction with the field can- 
not be treated perturbatively. Rather one must make the 
resonant approximation by finding a solution to the BdG 
equation (||) that is a mixture of the two adiabatic states, 
^a(r) and ^(r), in resonance: 



(3) 



*a0(r,t) = exp -i/2h / [E a (t) + E (t)]dt 



i(u/2)t 



b 2 s Vp(r)e 



-i{u/2)t 



Here the index s — 1, 2 labels independent solutions cor- 
responding to different initial conditions, b r s (t = 0) = 
6 Sir . After substituting the wave function (j|) into the 
BdG equation (^|) and averaging over the fast oscilla- 
tions, we readily find equations which describe the time 
evolution of the coefficients b 1 ^ 2 . Introducing a vector 
coefficient h s (t) we get 



ib s (t) 



5uj a p(t) 



V r 



a/3 



b.(t), (4) 



(2) where Sco a p(t) = us a p(t) —lu^lu and er + = (a x + ia y )/2, 



a x . y being Pauli spin matrices. A problem similar to 
Eq. (^) was first discussed by Landau and Zener in con- 
nection with molecular pre-dissociation . It follows 
from their theory that the asymptotic solution of (^) valid 
at t > i res is \b T s \ 2 = W a p + S rs (1 - 2W a/3 ) . The proba- 
bility W a f3 for the system to be in a different state after 
passing through the resonance, 



W a a = l- e 



(5) 



is determined by the product of the interlevel transition 
frequency \V a p\/h and the characteristic time St vcs = 
l^a/j^rcs)! -1 / 2 the system spends on resonance, 



7a/3 = \V a p/h\ \u af 3(t Tes )\ 1/2 . 



(6) 



If the characteristic time is long, 7 ^> 1 (St ICS ^> H/V a p), 
we have \b r s \ 2 ~ (1 — S r , s ) and the transition between 
modes a and (3 is complete. 

In the absence of relaxation, the Josephson current 
is solely determined by the dynamic evolution of the 
electron- hole wave function \I/(r,i) discussed above. Be- 
fore the high frequency field is switched on, the Andreev 
levels are in equilibrium and the Josephson current has 
its equilibrium value Iq. The deviation of the time aver- 
aged current, SI = I — Iq, is found using the asymptotic 
result for the coefficients b r s (t): 



51 = 



E 



W a fi (n a - Up) I (Ip - I a ) (2 - 5 Ea ,-Ep) 



(7) 

where I aj p = {2e/Ti)dE a> p/d(j>, and n a .p is the Fermi 
function, n Qj( g = np(E a ^). The summation in (^) is 
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over pairs of levels, {a, /3} re s, which have passed through 
a resonance ||. As we can see from Eq. (]?]) only transi- 
tions between states with energies of opposite signs (the 
Fermi level is at zero) contribute to SI. If the Joseph- 
son current through the channel is carried by a single 
pair of Andreev levels, the direction of the Josephson 
current is reversed for W > 1/2. It is essential to note 
that for t 3> ires there is no resonant coupling with the 
electromagnetic field to maintain the inverse population 
of Andreev levels responsible for the reversal of the cur- 
rent. Any inelastic relaxation mechanism will bring the 
system back to its equilibrium state and restore the orig- 
inal direction of Josephson current. This process of dy- 
namically switching on the inverse population of Andreev 
states and the subsequent inelastic relaxation manifests 
itself as a 'Somersault' of the Josephson current. 

Equation (0) is significantly simplified if the micro- 
constriction joining the two superconductors is an adia- 
batically smooth ballistic channel with quantized energy 
levels (modes) corresponding to the transverse electron 
motion. In this situation each propagating mode pro- 
duces one pair of Andreev levels of the type (P — one 
level above- and one below the Fermi energy — and the 
transmission eigenvalue T n reduces to the normal elec- 
tron transmission coefficient of the n-th mode, 



T„(fcf 



AR 



(1-RY 



■ sin 2 (k n L) 



(8) 



where k n = \/k 2 F — w 2 n 2 /d 2 , and R is the probability for 
electron backscattering at the Sm-S boundary. The elec- 
tromagnetic field couples levels only within these pairs, 
i.e. if a = ± one has V aa ' = &n,ri&a.-<j l Vn and conse- 
quently uj aa i = 8 n _ n '5 a ^ a iLj n . The matrix element V n 
has the form: 



\V n 



T 2 

V 2 —T 3 



R cos 2 (k n L) sin 2 </> 
£sV"(l-^) 2 l-T„sin 2 (0/2)- 



(9) 



In the ballistic case the total current at zero temperature, 
taking into account Eq. (f^), can be written in the form: 



I = ^I n (l-2W„), 



(10) 



where the summation is formally done over all transverse 
modes and W n is assumed to be zero for those Andreev 
levels which have not passed through a resonance. 

As one can see from Eq. (||), the electromagnetic field 
can only mix Andreev states if the probability for elec- 
tron backscattering at the Sm-S boundaries is finite. A 
'Somersault' which changes the total momentum of the 
system can evidently only occur if backscattering is possi- 
ble. It is interesting that if there is a geometric resonance, 
k n L = 7i7T, the Somersault effect still takes place ac- 
cording to (j^) even though the microconstriction is fully 
transparent for electrons in mode n, T n (kp) = 1- This 



fact provides an important possibility of distinguishing 
between resonant transmission through the double bar- 
rier structure and balllistic propagation in the absence of 
normal backscattering from Sm-S boundaries. 

It is also interesting to consider the case when there 
is no drift with time of the Andreev levels, i.e. when 
Vq(£) = 0. A permanent resonant coupling of Andreev 
states is now possible and leads to an equalization of the 
average level population [jfl). Solving equations (||) we 
find a result for the Josephson current that deviates from 
its stationary value: 



5> 



\8uj r , 



(2v n /ny 



(ii) 



In particular we find that resonant coupling, i.e. when 
Sui n = 0, results in a complete blockade of the Josephson 
current. This phenomenon can be observed as a sequence 
of dips if the Josephson current is plotted as a function of 
the frequency of the electromagnetic field. The number 
of dips should equal the number of propagating modes 
in the ballistic constriction and the dip amplitude should 
equal the single mode contribution /„ to the stationary 
Josephson current. 

In the discussion so far we have assumed that the phase 
difference between the two superconductors is constant in 
time. For the Somersault effect to develop as described 
above it is necessary that </> be fixed to a high degree of ac- 
curacy, otherwise the system will escape from resonance. 
This fact implies that the change of the interlevel dis- 
tance due to the phase shift 8(f) should be smaller than 
the nonlinear resonance width: %(du) a p/ d(j))5(j) <C V a p. 
To estimate 6<j) we first observe that <fi is related to the 
shielding currents induced by the enclosed magnetic flux 
in the SQUID geometry. This Meissner current flows 
near the surfaces in the thick part of the ring shown in 
Fig. 1. The diameter of the ring cross section is d 3> A, 
where A is the penetration length of the magnetic field. 
The current flowing through the small ring segment that 
forms a narrow channel is much smaller than the Meiss- 
ner current in the main part of the ring. The Somersault 
effect, which affects this small part of the induced cur- 
rent, will result in a redistribution of Meissner current 
in the bulk ring region. The corresponding phase shift 
will be of order S<fi ~ (£/£s)(l/JVx), where C is the ring 
length, £5 is the superconductor coherence length, and 
N\ ~ k 2 Xd is the number of transverse modes carrying 
the Meissner current (here ko is the Fermi wave vector 
in the superconducting ring) . The criterion for the phase 
difference to be stable enough to observe the resonance 
is thus (V a p/A)(L/£) » (1/JVa). 

When the Josephson junction is connected to a cur- 
rent source, the current rather than the phase differ- 
ence is held constant. In this situation another inter- 
esting phenomenon occurs; the Somersault effect causes 
phase slips at times when Andreev levels are resonantly 
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coupled. The phase slips in turn manifest themselves 
as voltage spikes. The characteristic time scale of a 
phase slip is determined by the Josephson plasma fre- 
quency wj = [(A/^)(l/i? C J )] 1 / 2 , which for the ge- 
ometry of Fig. 1 and Rq = Ti/e 2 can be written as 
luj = (A/h)(£sL/ 'd 2 ) 1 / 2 . This relationship allows us to 
estimate the amplitude and the duration of the voltage 
spike as U — (Ti/2e)(/> ~ {%/e)u>j and At ~ 1/ojj respec- 
tively in the case where the duration of the interlevel 
transition, St rcs , is shorter than At i.e. if l/u)j 3> St Tes . 
In the opposite case, 1/toj <C St YOS , a more complex 
analysis of the joint dynamics of the phase- and inter- 
level transition gives the estimates U ~ V^p / eA 1 / 3 , and 

At~h/V^ 3 A^. 

The Somersault effect discussed above relies on the 
preservation of phase coherence during the dynamical 
evolution of the Andreev states and the redistribution 
of level populations at resonance. This implies that the 
phase breaking time t associated with inelastic relaxation 
should be the largest time scale in the problem: 

l/w a fi ~ h/A < %/V a p < St ICS < t, 1/uj < r. (12) 

Since the main part of the Andreev state wave function 
is localized in the superconducting electrodes the most 
important relaxation mechanisms are recombination pro- 
cesses due to electron-phonon interaction in the super- 
conductors. According to Ref. |l2| a typical estimate of 
the recombination time in a superconductor is At/Ti > 
10 2 . This value makes it possible to fullfil the inequalities 
( |l2| ) in a device with normal region length L w £s m if the 
amplitude V u is of order of 10 /iV. Meanwhile, an ampli- 
tude of the voltage spikes in a current biased junction is 
of order of 1-10 /iV. 

In conclusion, we have shown that a microwave elec- 
tromagnetic field can drastically change the Josephson 
current carried by quantized modes in a mesoscopic con- 
striction, giving rise to sequential reversals of the current 
in a phase biased junction and to voltage spikes in a 
current biased junction. Observation of the effect, which 
seems to be possible in currently available gated semicon- 
ductor heterostructure devices, would be a direct mani- 
festation of the discrete nature of the Josephson current 
in a quantum constriction. 
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FIG. 1. Sketch of a gated semiconductor heterostructure 
in which a two-dimensional normal constriction is confined 
between superconducting electrodes. Inset: The phase differ- 
ence between the two superconductors can be kept constant 
in the SQUID configuration shown (see text). 

FIG. 2. Microwave induced transitions and currents car- 
ried by the Andreev levels in the normal constriction of Fig. 1 
during a slow drift of the levels through resonance: (a) before 
reaching the resonance; (b) during the resonance; (c) after 
passing the resonance. Arrows show the direction of current, 
dashed arrow correspond to unpopulated level. 
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